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A semi-classical, many-body atomic model incorporating a momentum-dependent 
Heisenberg core to stabilize atomic electrons is used to study antiproton capture on 
Helium. Details of the antiproton collisions leading to eventual capture are presented, 
including the energy and angular momentum states of incident antiprotons which 
result in capture via single or double electron ionization, i.e. into [He ++ p or He + p], 
and the distribution of energy and angular momentum states following the Auger 
cascade. These final states are discussed in light of recently reported, anomalously 
long-lived antiproton states observed in liquid He. 
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I. INTRODUCTION 



The "trapping" of antiprotons stopped in liquid He into long-lived states has been re- 
ported recently by Iwasaki et al. Approximately 3.6% of the incident antiprotons exhibit 
delayed decay, of the order of microseconds compared with picoseconds for prompt decays, 
opening possibilities for further experimentation, including antihydrogen production. 

The initial slowing and capture of in-flight antiprotons is via energy loss due to Auger 
excitation and ionization of the atomic electrons, while the later stages are dominated by 
radiation, and eventual annihilation. Auger transitions are characteristically many orders 
of magnitude more rapid than radiative decay, making the dynamics of the Auger process 
central to understanding the delayed antiproton annihilation. 

Yamazaki et al. @ , following the work of Condo || and Russell ||, proposed that 
the long-lived antiprotons were captured into "meta-stable, circular states," of principal 
quantum number n = y M*/m e m 38 and with i w no — 1, of the exotic neutral He + p 
atom. (M* = |Mp is the antiproton reduced mass with respect to He and m e is the electron 
mass.) 

From these states further Auger transitions would be highly inhibited: 

1) A£ > 1 is highly unfavored and the atomic excitation energies are large compared 
with the antiproton spacing for An = A£ = — 1, and 

2) The Stark effect, which normally admixes £ = states due to the presence of other 
atoms in the liquid and enhances p annihilation in the nucleus for atoms stripped of inner 
electrons, is suppressed due to the removal of ^-degeneracy by the presence of the remaining 
electron. 

Thus decay is limited to the much slower radiative transitions, which have small photon 
energy, and antiproton annihilation is delayed. 

In order to study the capture of antiprotons by He atoms, we modeled antiproton colli- 
sions with He using a classical description of the antiprotons and a semiclassical description 
of the electrons, in which the electrons of the target system are stabilized by a momentum- 
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dependent Heisenberg core [5]; the Pauli principle can be ignored here since the atomic 
electrons are in antiparallel states. 



II. QUANTUM AND ADIABATIC SOLUTIONS OF THE HE-p SYSTEM 

Consider now the quantum states of this system, an exotic "molecule" consisting of the 
He nucleus, i.e. an a particle, the antiproton and N =0, 1, or 2 electrons. Working in 
atomic units, with h = m e = e = 1 (ignoring the negligible effect of electron reduced mass), 
the Hamiltonian is given by 
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where r*j and R are the electron and antiproton coordinates relative to the helium nucleus, 
and pi and P are the corresponding momenta; M* = 1469 is the reduced mass of the 
antiproton in atomic units. 

The zero electron (double ionization) problem is simply the He + ion with the p replacing 
the electron. The energy levels are given by 

Z 2 M* 2M* 



Note that n = a/M* ~ 38 corresponds to a p orbit of the same size and energy as the 
Is electron orbit; because of the large M* and consequent large quantum numbers for the 
antiproton states of interest, the antiproton may be treated classically. 

The one electron (single ionization) problem is described well in the adiabatic Born- 
Oppenheimer approximation. We note that this two-center, one-electron problem is sepa- 
rable in prolate spheroidal coordinates and can thus be calculated to high accuracy, cf |3J, 
||. In terms of the coordinates defined above, the energy levels e n (R) of one electron in the 
potential of the He nucleus and fixed antiproton are given by 
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and are the potential energy for the Born-Oppenheimer antiproton eigenvalue equation, 
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Here P 2 = —R~ 1 (d 2 /dR 2 )R + L 2 / R 2 includes vibrational and rotational excitation. These 
Born-Oppenheimer states are not stationary states of the full problem, however, since there 
remain velocity terms in the full Schroedinger equation which couple the adiabatic levels. 
Quite generally we can expand the total wave function in terms of the eigenfunctions of 
Eq. (§. 

^(VisionRes. , R) = ^ $ n {R)ipn(VisionRes. , R) . (5) 
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The Schroedinger equation for the $ n is then 
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where \n> corresponds to ip n and the bra-ket integration is over VisionRes. . 

In a time-dependent calculation, these coupling terms lead to Auger transitions between 
the adiabatic states. Note that the first coupling term is essentially of the dipole form (the 
operator is P). Thus transitions with \A£\ > 1 are (progressively more) inhibited. 

Alternately, Yamazaki and Ohtsuki have obtained approximate solutions to the Born- 
Oppenheimer problem using configuration-mixing techniques, in which antiproton states are 
calculated with the electron in the ground state, and then mixed with excited states of the 
electron to obtain the system levels shown in Fig. 1. Here, the system levels L = N — 1 are 
proposed as the boundary of the allowed E, L states, with the states with L > 31 proposed 
as the region of metastability. 

Because of the infinite range of the central Coulomb potential, the Born- Oppenheimer 
energy level diagram of E(L), where L is the total angular momentum, in theory extends to 
the right indefinitely; this can be seen by considering first a simple ap state in the energy 
region covered in Fig. 1. An electron can then be placed in a hydrogenic orbit of arbitrarily 
small energy and arbitrarily large angular momentum without disturbing the p. Only some 
of these states were observed in our calculations, however; electrons excited to higher angular 



momentum states during the collision process also gained enough kinetic energy to escape 
from the (exotic) atomic system. 

Capture of p into the two electron negative ion was also not observed in our calculations; 
this negatively charged system appeared only briefly during the initial stages of antiproton 
capture, and was unstable to electron ionization during the antiproton capture and decay 
process. 



III. THE SEMI-CLASSICAL MODEL 

In the semi-classical model ||, the Hamiltonian of the undisturbed Helium atom is 

2 2 n r 2 ri 2 

with Vh the momentum-dependent Heisenberg core which prevents the collapse of the elec- 
trons into the He nucleus, given by 

V H = ^ r2 exp{a[l-(r P /^}} (8) 

We choose || a — 1.0, =2.767. The Hamiltonian then minimizes at an energy 
of —2.78 with r\ = r 2 = 0.63, in fair agreement with the exact ground state energy of 
—2.9037 • • ■ at an electron mean radius of < r >^ 0.59. 

The total Hamiltonian is then 
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To model collisions, Hamilton's classical equations of motion 

dxi dH sc dpi dH sc 



+ V H (ri,pi) 
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dt dpi dt dxi 

are solved for VisionRes. i, pi, VisionRes. 2, P2, R, P , the coordinates describing our 
semi-classical system, and integrated over time using the venerable ordinary differential 
equation routine ODE. 
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IV. COLLISION CALCULATIONS 



Monte Carlo calculations of antiproton collisions with He were performed as follows: 

(1) The target He atom was prepared in its ground state with random orientation and 
parity inversion of the electron coordinates and momenta, as described in 

(2) The projectile p was launched at this target atom with an initial energy of about 80 
eV (about 3.0 a.u.) and with an impact parameter b randomized with equal areas (ndb 2 ) up 
to an energy-dependent b max . 

(3) Sequences of collisions were followed from one encounter to the next until capture 
occurred. 

Our initial energy and b max were chosen so that collisions would never result in antiproton 
capture on the first encounter, and the final energy of the antiproton after each collision 
could be used as the starting energy for a subsequent collision, again with Monte Carlo 
initial conditions for the ground state of the new target and for the new impact parameter. 

As the antiproton slowed, b max was increased in stages: we began with b max = 1.0; when 
p energies dropped below 2.3 a.u., the level at which small numbers of captures first began 
to be observed, b max was increased to 2.0; when p energies dropped below 1.2, somewhat 
above the range where capture cross sections began to flatten out with b max of 2.0, b max was 
increased to 3.0. 

While somewhat artificial, by selecting for collisions during which something happened 
this method allowed us to study a realistic distribution of antiproton projectiles during the 
last stages of antiproton cascade down to eventual capture, typically after five to twenty 
sequential collisions, while decreasing the amount of computer time spent modeling weakly 
interacting collisions of large energy and impact parameter. 

A total of 4,000 collision sequences (approximately 41000 individual He- antiproton col- 
lisions) were followed. 916, or about 23%, ended with capture via single electron ionization 
into the neutral exotic aep, while the remaining 77% resulted in capture via double ionization 
into a positively charged ap 
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Fig. 2a shows our total energy loss cross sections, d(aAE)/dE, of the p as a function 
of initial p energy. In the higher energy collisions at the start of a collision sequence, 
little slowing of the antiprotons occurs for large impact parameters. As the antiprotons 
slow, energy loss increases at larger impact parameters until the antiprotons drop into the 
energy range from which capture begins to occur, slightly below the energy at which we first 
increased the maximum impact parameter; step increases in these energy loss cross sections 
are clearly visible at energies of 2.3 a.u. and 1.2 a.u., where we increased the maximum 
impact parameter. 

Fig. 2b shows our antiproton capture cross sections for these collisions via single and 
double ionization, again as a function of antiproton energy. As the antiprotons slow, the 
maximum impact parameter at which capture occurs grows; since the central Coulomb 
potential of the He nucleus has the form 1/r, the total capture cross-section becomes infinite 
for very low energies. Here, our total capture cross section levels out to Ttb 2 max at the lowest 
antiproton energies, where all the antiprotons are captured; our maximum impact parameter 
of 3.0 for these lower energy collisions was chosen so that this saturation occurred only for 
p energies of approximately 0.1 a.u. or lower. 

Due to the limited amount of energy and angular momentum that can be exchanged by 
the massive antiproton and the bound electrons of the He atom, antiproton energies and 
angular momenta which result in capture fall into a fairly narrow range, as shown in Figs. 
3a and 3b. Fig. 3a shows the range of incident antiproton energies and angular momenta 
which result in capture via double electron ionization; the adjacent range of higher energy 
and angular momenta antiprotons which were captured via single ionization is shown in Fig. 
3b. 

Fig. 4a shows the energy and angular momentum of the final states into which the an- 
tiproton was captured via double ionization. For capture via single ionization, Fig. 4b shows 
the ape system states; in Fig. 4c, just the antiproton states of these systems are plotted. 
Here, after capture the antiproton continues to interact and exchange energy and angular 
momentum with the remaining electron, so each of states plotted in Fig. 4c represents a 
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time average of the antiproton configurations. Fig. 5a,b, showing the changing electron and 
antiproton radii as they interact over time, is an example of a typical aep system dynamics 
after antiproton capture. 

Density contours of these scatter plots, in percent of total collisions per unit E per unit 
L, are combined in Fig. 6, illustrating how the adjacent bands of incident antiprotons of 
Fig. 3 drop into the captured states of Fig. 4, and the relationship of the states of the 
antiprotons captured via single ionization to the states of the total system. 

Also of interest are the transition E and L distributions of the antiprotons when they are 
first captured by the He atom. Fig. 7a shows the initial capture states of the antiprotons 
in what will become the doubly ionized ap, when the system was an unstable aep or, in 
some cases, aeep. Fig. 7b shows the higher angular momentum states of the antiprotons 
captured into what will become the singly ionized aep, again at the point of initial antiproton 
capture. Comparison with Figs. 3 and 4 shows the incident antiprotons with higher angular 
momentum of Fig. 3b clearly separating into the higher angular momentum states associated 
with the more stable aep of Fig. 4b, relative to the lower angular momentum antiprotons 
in those systems that will quickly ionize to the ap states of Fig. 4a. 

V. CAPTURE ANALYSIS 

The quantum antiproton levels for the simple ap are given above by Eq. 

Since no ~ 38 has generally been considered the lower limit of the Auger cascade, i.e. 
the level below which p orbits are progressively less disturbed by any electrons which may 
still be attached, this has been assumed to be the boundary between Auger and radiative 
decay, where exotic projectiles are expected to accumulate following the Auger cascade 
before further decay by radiative processes; much below this level decay proceeds primarily 
by radiative capture irrespective of the state of ionization. 

In our calculations, antiprotons captured into ap are distributed in this region; Fig. 4a 
shows the antiprotons stacked up against the classical L = n centrifugal boundary in the 
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region of n « 30 — 40, with lower L states more heavily populated in the region n < 35. It 
is not surprising that low energy antiprotons are captured into states of n < n via double 
ionization, since the antiproton must exchange enough energy with the atomic electrons to 
boost both of them into (positive energy) escape from the He nucleus. 

It is also generally assumed that the final populations of these states are distributed 
according to the statistical factor 2£+ 1. Fig. 8a shows the distribution of the ap states vs. 
L/n, where the classical circular angular momentum is given by 



In this figure, the total L-state distribution, averaged over the different n-states, is approx- 
imately linear in L. 

The separate higher energy and angular momentum band of aep system states shown in 
Fig. 4b can be compared with the aep system level diagrams proposed by Yamazaki et al, 



additional angular momentum contributed by the remaining electron in the singly ionized 
systems are indeed populated in our calculations. 

The distinct, narrow band of high L excited p states present in the singly ionized system 
(Fig. 4c) stack up against a modified centrifugal barrier, at a lower value of L than in the 
doubly ionized system due to the screening effect of the remaining electron. The degree to 
which the antiprotons stack up against this barrier is illustrated by the population distri- 
bution of Fig. 8b. Here, the population distribution is plotted against L/n e ff, where ra e //, 
labeling the energy levels of the screened system, 
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0, shown in Fig. 1; as discussed above, some of the higher L system states resulting from 
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is given by 
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with Z e /f, defined in terms of the central radial force felt by the antiproton 



F r = ~ 




(14) 
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VI. METASTABILITY 



As discussed above, it is a small percentage of antiprotons in stable, high L "circular" 
states, from which the antiproton cannot decay to the point at which it overlaps the nucleus 
and annihilation occurs, that has been proposed as the source of the observed metastability 
in long-lived exotic aep atom In our calculations, a relatively large number of the aep 
atoms are in these higher L states; it is, in fact, the higher initial angular momentum of 
the incident antiprotons which results in their capture via single ionization, as can be seen 
by comparing the initial antiproton angular momenta, at the start of the collisions resulting 
in antiproton capture, of Figs. 3a (double ionization) and 3b (single) with the final system 
angular momentum distributions shown in Figs. 4a and b. 

The question of what constitutes a stable, singly-ionized aep system is of course related 
to how long the system is followed. In our initial calculations, systems were followed for a 
maximum time of 25,000 a.u., primarily due to limits in computational resources. Indeed, 
some small fraction of the systems classified as singly ionized have very small binding energy, 
as can be seen in Fig. 4c. In most of these cases where the p is very loosely bound, the p will 
eventually escape from its exotic system, to be captured by the next He atom it encounters. 

To further investigate aep stability in our model, 500 systems were followed to a maxi- 
mum time of 10 5 a.u. As shown in Fig. 9, the percentage of singly ionized systems surviving 
on these time scales is well- described by a double exponential of the form 

Paep(t) =Po + Pie' t/Tl + P 2 e- t/T2 = [22.3 + 82.6e-' /2051 + 7.67e~ t/57 ' 650 ] % (15) 

This fit gives a survival of 22.3% up to time 10 5 a.u. compared with the observed metastable 
population of 3.6% for survival beyond about 1/xs ~ 4 10 10 a.u. Some decay in the singly 
ionized population is observed at the end of this time scale, suggesting the existence of 
additional, longer time scale decay constants for this population. 
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VII. CONCLUSIONS AND FUTURE 



Semi-classical Monte Carlo calculations of the collisions of antiprotons with He offer new 
details of the slowing and capture processes, including the dependence of the final products 
on the incoming energy and angular momentum and the state distributions of the singly 
and doubly ionized exotic aep and ap at the end of the Auger cascade. These details shed 
light on the recently reported metastable states of aep. 

Further questions to be investigated with this method include the effects of impurities 
in the target medium on the slowing and capture process, and more detailed study of the 
Auger process from which the high L meta-stable states arise. The population decay of 
these systems near the end of our time scale invites further study of the dynamics of this 
decay, and of their relation to the much longer time scales of the metastability discussed 
by Yamazaki et al, ||, who quote mean lifetimes of about 3 /is ~ 1.25 10 11 a.u., which is 
beyond the range of our current computational tools. 

We can suggest some mechanisms which might reduce survival to 1/zs (which we note 
is 4 10 5 times our longest runs. One is Stark mixing into s-states due to the other helium 
atoms in the liquid, leading to annihilation on the He nucleus. This is inhibited by the one 
electron, but may still be significant. Another is annihilation from the higher L-states which 
we find to survive. We have the population of these states. These could be used by those 
who have calculated the Born-Oppenheimer states to calculate the annihilation rate from 
their p wave functions at the He nucleus. 
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FIG. 1. The energy levels of of the aep system courtesy of Yamazaki and Ohtsuki ||. 

FIG. 2 a. Energy loss cross sections for 41,417 p-He collisions, as a function of initial p 
kinetic energy, T. 

FIG. 2 b. Capture cross sections for 4,000 p-He collisions, as a function of initial p kinetic 
energy, T. Capture into pea states are deno ted by (+), into pa by (x) and total by (*). 

FIG. 3a. Initial state of the p for the last pass in the cascade before capture into a, pa 
state. 

FIG. 3b. Initial state of the p for the last pass in the cascade before capture into a pea 
state. 

FIG. 4 a. Final energy and angular momentum for 3,084 ap states. 

FIG. 4b. Final total system energy and angular momentum for 916 aep states. 

FIG. 4c. Time averaged p energy and angular momentum for 916 aep states. 

FIG. 5 a. Radial coordinate of the p as a function of time for a particular run, T = 
0.5, b = 0.977. 

FIG. 5 b. Radial coordinate of the electron as a function of time, in comparison to a 
portion of the antiproton trajectory shown in Fig. 5 a. Because of the rapid oscillation of 
the electron, the figure does not show all of the calculated points. 

FIG. 6 Density contours, in percent per unit E per unit L. The various regions (islands) 
describe: A. The state of the incident antiprotons for the final collision which ends in capture 
via double ionization. B. The state of the incident antiprotons for the final collision which 
ends in capture via single ionization. C. Final (time-averaged) p E and L in the aep 
configuration D. Final system states for the ap configuration, and E. Final system states 
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for the aep configuration. 

FIG. 7 a. Energy and angular momentum at initial p capture for collisions leading to ap 
products. 

FIG. 7 b. Energy and angular momentum at initial p capture for collisions leading to 
aep products. 

FIG. 8 a. Distribution in the ratio of angular momentum states L to the classical max- 
imum angular momentum state L max = n for final ap states. Averaging over the different 
n-states shows an approximately linear distribution in the L-states populated. 

FIG. 8 b. Distribution in the ratio of p angular momentum states L to the classical 
maximum angular momentum state L max = n e ff for final aep states. Averaging over the 
different n-states here shows the concentration in high L states of the ps trapped into aep 
systems. 

FIG. 9. Survival as a function time of aep systems after initial p capture. 
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